module probabilidades
 use parametros
 contains
 function prob(ket0,j,k) 
  real*8::prob
  complex*16, dimension(1:4,-dm:dm,-dm:dm)::ket0
  integer::j,k
  prob=abs(ket0(1,j,k))**2+abs(ket0(2,j,k))**2+abs(ket0(3,j,k))**2+abs(ket0(4,j,k))**2
 end function prob
end module probabilidades
